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Abstract 

In response to the experimental discovery of double perovskite Ru- base material, 
The magnetic and electronic properties of the double exchange mechanism are studied 
theoretically. The electron band structure and the spin wave spectrum response at 
zero temperature with the canting background are computed from the first principle 
Hamiltonian. The long range Coulomb interaction between charge carriers is found 
to be crucial in stabilizing the canting order. A Monte Carlo simulation is performed 
on de-Gennes's effective magnetic Hamiltonian with classical spin approximation. The 
peaks of the specific heat as a function of temperature, which correspond to FM(AF) 
and canting transitions, and its response to the external magnetic field are consistent 
with the observation. For certain region of parameters, a more complex canting order 
than layer type is found at sufficiently low temperature. The effective interaction 
between carriers via exchanging a magnon wave is found to be rather similar to that 
via exchanging an acoustic phonon. The pros and cons to a genuine superconducting 
order in this system is analyzed. 



1. Introduction 



Recently, a new type of materials of double perovskite structure is syn- 
thesized by [1], which displays a number of interesting physics [1-2]. The 



1 



molecular formula are Sr2Y Rui-xCuxOq, Ba2Y Rui-xCuxOq, etc.. The 
magnetic properties are unusual with two phase transitions as a function of 
temperature. The underlying mechanism was conjectured to be the double 
exchange effect [3-5], the same mechanism partially responsible to the colos- 
sal magnetoresistance in the single perovskite compound, La2-xC axMnO 4^ 
[6]. In addition, the materials become superconducting for x > 0.04 and its 
transition temperature reaches Tq = 20 — 30K for x = 0.1. The coherence 
length is about ^ = 35A and the super phase possess several novel features 
different from cuprates. 

The double perovskite crystal structure of the parent compound {Sr, Ba)2YRuOQ^ 
can be approximated by a cubic structure with the side length of a funda- 
mental cube / ~ 8.16A. The {Ru,Cu) ions are located at the face-center 
sites, the Y ions are sitting at the octahedral sites and Sr ions at tetrahedral 
sites. The O ions arc placed at the middle of the shortest Y — Ru bond. A 
small fraction of Ru is replaced by Cu through doping. 

The charge balance of the ions Sr{Ba), Y, Ru and O of the parent com- 
pound are +2, +3, +5 and —2. Except Ru^^, the outer electronic configura- 
tions of all other ions are closed shells of the continuous rotation group. The 
configuration of Ru^^ is AS. It is likely that the strong Hund's rule coupling 
enforces the spins of the three d electrons parallel with total spin S* = 3/2 
and their orbital wave function fill up a multiplet of the T2 representation 
of the cubic group, supported by xy, yz and zx (notably, the lobes of these 
wave functions point to the directions of the twelve nearest neighbors of a 
face center site). With only the nearest neighbor antiferromagnetic exchange, 
the Neel state may assume a layer structure [7] perpendicular to one of the 
4-fold axis with the spins of Ru~^^ alternating along successive layers. This is 
consistent with the tetragonal symmetry observed by X-ray diffraction data 
[1]. Replacing part of Ru ions offset the charge balance and creates itinerant 
holes with density 2x per molecule of {Sr, Ba)2YRui-xCuxOQ. This may 
also break the symmetry further to an orthohmobic structure as indicated 
by the X-ray data of the doped compounds. 

With itinerant holes, the state of a Ru ion is a superposition of states 
with charge -|-5, as in the parent compoimd and states of charge +6. In the 
latter case, Hund's implies a total spin S = 1. Regarding the latter states as 
combinations of the product states of an ion core Ru^^ of spin S — ^ and a 
hole of spin |. Houd's rule will prefer antiparallel spins of the ion core and 
the hole. The Hamiltonian combining both the antiferromagnetic exchange 
and the double exchange reads 

H = J2 JiA-Sj -\Y. + + ^ S^^]o^j + UcaaU (1-1) 

ij ij j 
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where 5^^ is the total spin operator of the ion core, are the creation 

and the annihilation operators of the itinerant holes, each of two components 
with spin up and spin down, Jy denotes the antiferromagnetic exchange, 
tij denotes the hopping amplitude of the itinerant electrons and A denotes 
the exchange interaction of the itinerant electrons with the ion cores. The 
Coulomb interaction, Ucoul.i of (1-1) is given by 

(1.2) 

where the first term represents the on-site Coulomb repulsion, e denotes the 
dielectric constant of the ion cores and :(...): enforce the normal ordering. 
As we shall see, the long range Coulomb interaction is rather important 
to maintain a homogeneous canting magnetic order against clustering. The 
magnetic anisotropicity, the possible Jahn- Teller efi'ect due to the degeneracy 
of d-band, as well as the impurity effect, which might be relevant in real 
systems, have been neglected for simplicity. 

It was pointed out by Anderson and Hascgawa [4] that the double ex- 
change energy, steming from intra-atomic Hund's rule coupling, ought to be 
much stronger than the hopping energy, i.e., \S » t. In these circumstance, 
the effective magnetic Hamiltonian after integrating out the electronic de- 
grees of freedom can be approximated by [6] 

ij ij 

with 



for low carrier densities 

This paper addresses mainly the theoretical issues of the magnetic prop- 
erties associated with the double exchange effect. Throughout this paper, 
we shall adopt tetragonal order for the magnetic states, which preserves the 
maximum subgroup symmetry of the original cubic group. Only J^j's and 
tjj's for nearest neighboring bond are kept different from zero for simplicity 
and we introduce the notations that 

J _ f interlayer bond < ij >; (15) 

1 J', for intralayer bond < ij > 
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and 



^ _ ( t for interlayer bond < ij >; (1 Q) 

\ t', for intralayer bond < ij > 



Accordingly, the parameters bij for the effective magnetic Hamiltonian (1.3) 
read 

b for interlayer bond < ij >; 



1 6', for intralayer bond < ij > ' v ■ / 

On a cubic lattice of A/" sites with a periodic boundary condition, there are 
2J\f intra-layer bonds and 4J\f interlayer bond. Following de Gennes [5], we 
a homogeneous spin orientation in each layer and a canting order between 
successive layers. The canting angle as follows the minimization of (1.3) is 

This paper is organized as follows. In the section 2, we shall expand the 
Hamiltonian (1.1) around a canting order and examine the band structure of 
the itinerant holes. The first principle calculation of the magnon spectrum 
at zero temperature will be displayed in the section 3. The Monte Carlo 
simulation of the magnetic Hamiltonian (1.3) at arbitrary temperatures is 
presented in the section 4. The section 5 is devoted to a sketch of possible 
coupling between the magnon and a superconducting order. In the section 
5, we shall come back to the real systems and discuss the possibility of a 
genuine superconducting order in them. 



2. Mean-Field Expansion of the Hamiltonian around Canting 



Denote the side length of a basic cube by I and the cubic axes, (100), 
(010) and (001) by x,y,z. The locations magnetic ions are given by 

-R = mei + 7126*2 + 7^363, (2.1) 
where ni, n2 and 713 are integers, 

ei = ^(y + 5), (2.2) 

e-2 = l(5 + x), (2.3) 

e3 = ^(^ + y), (2.4) 



(a) (b) 
Figure 1 (a) The magnetic ordering. The magnetic sites of the sublattices A and B are 
marked expUcitly. (b) The spin frame. The angle between (a, and Cs axes is O. 

and the total number of sites is J\f. Adapting the tetragonal layer magnetic 
ordering, we divide the lattice into two sublattices shown in Fig 1. with 
ni + n2 = even on the sublattice A and ni + n2 = odd on the sublattice 
B. The expectation value of the ion spins of the sublattice A is denoted by 
Sa — S(a and that of the sublattice B hy Sb — S(b- We choose the C-axis 
along the direction of Sa + Sb and the r^-axis as 

— * — * 

Ca X Cb ,^ c-x 
sm fc) 

with the mutual angle between Sa and Sb- We introduce further 

I = ^ X C, (2.6) 

= r) X U (2.7) 

and 

Is = X Cs. (2.8) 

The sets (^,?7, C); (^^;^;C^) and {^BjVXb) each forms a right hand coordi- 
nate system. The first of them will be referred to as the spin frame. Since the 
magnetic anisotropic energy has been neglected, the orientation of the spin 
frames does not couple with the crystal axis. The relevance of the magnetic 
anisotropy will be discussed in the Appendix A. 
With large spin approximation [8] , we write 

Sj = V f (cj + c])iA - iJ f (c,- - 4)77 +{S- c]cj)CA + ... (2.9) 
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for the site j on the sublattice A and 



2 {dj + d])iB - -{dj - +{s- didjXB + 



(2.10) 



'-'3 ~ 

for the site j on the sublattice B, where cj, Cj, dj and dj satisfy the commu- 
tation relations of elementary boson operators. After some algebra, which is 
shown in detail in Appendix B, the Hamiltonian (1.1) is expanded as 



where 



H — Hq + Hi + H2 + + H4, 

Ho = 27yA52(j' + 2Jcose), 
Hi = 4JS^V2AfiCo + Co) sin 6, 



(2.11) 

(2.12) 
(2.13) 



e 



- J S siv? ^ vp{zpz_^ + zKzi - Cp<-p- - Clp<}) 
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with = 2S 

p 



p 

J'{-2 + u^) - 4Jcos 6 =F Jvp cos^ § 



+Zg4'\ ffPlT-ip^ q + -pir+V-,? + three-magnon terms, (2.15) 

P~i ^ P~i 



p,q 



.t./,t 



p+i. 



Ha^ 



A 



^ X f-^ p+5 



27V ^--^L -p'-i 



(2.14) 



+ ^ir_|V+f )'^^-+|/^i^3V'^_| + four-magnon terms (2.16) 



and 



+4°^'(pi)a/3(pi)a'/?'V'l£>t_. ,V',-,,g>V',^_£ 



(2.17) 



In (2.15)-(2.16), only the terms wc shall use are exhibited explicitly. The 
magnon operators and z^ of (2.13)-(2.16) are related to Cj and dj via 



(B.l), (B.2), (B.19) and (B.20) in the appendix. The hole operator i/j^ of 
(2. 14) -(2. 17) is a 4 X 1 column matrix and is related to the 2x1 column 
matrix \E'j through (B.3), (B.4) and (B.21). The 4x4 Dirac matrices 

n=f- lV^3=f- 1) (2.18) 

and 



with 

[Pa,n]=0 (2.20) 

and 

r± = ^{Ti±iT2) (2.21) 

have been introduced accordingly, where / is a 2 x 2 unit matrix and a^s 
are Pauli matrices. The Bloch momentum p(i/, q) throughout (2.13-17) cor- 
responds to the translational invariance within each sublattice. In terms of 
the dimensionless papameters 

ai=^^{Px+Py), (2.22) 

«2 = ^(-Px+Py) (2.23) 

and 

«3 = Pzl, (2.24) 

the Brillouin zone is specified as tt < a i, 0:2, 0:3 < n, the function and 
of (2.14) is given by 

= cos CKi + cos a2 (2.25) 

and 

ai Q!2 as 

v^=4cos — cos — cos — , 2.26) 
P 2 2 2 

and the bare Coulomb propagators D^^^ and D^S''' of (2.17) read 

40 / (2.27) 

9 i(2a2^ + Q;|) Z ^ ^ 

Df = C/-H!!i?. (2.28) 
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The Hamiltonian (2.11)-(2.17) is invariant under the tetragonal crystal 
rotation, inversion, and the C-transformation introduced in the appendix B, 

Cz^~^ = -zp (2.29) 

CCpC"' = Cp (2.30) 

and 

C^f-^ = -psipp. (2.31) 
The matrix of the hole kinetic energy is given by 

Ep- 11 = -t'up - tvpps (t3 cos ^ - n sin ^ j + ^AS'rs - fi, (2.32) 
which is easily diagonalized and produces four energy bands: 

= (2-33) 

for C = 1 and 

for C = —1, where s = ±1 and 



A± = J -X2S'2 ± A;5to„^cos - + M. (2.35) 
The wave functions corresponding to the four bands are 

^2.36) 

for ei„ and 

p,s 

^f,s=(^r) (2-37) 



p,s 



for e-.„, where 2x1 matrices d>^^ are defined by 



1 

A? 

V 



(^\S ± to^cos ^)cT3 T tw^i sin ^^%s = (2-38) 



The electron spin ceases to be a good quantum number and the only degen- 
eracies in Brillouin zone correspond the points related by the space inversion 



and the proper tetragonal rotations. The density of states of the four bands 
are plotted in Fig. 2. 

The free electron propagator with energy pQ and momentum p is given 

by: 

/oo 
dte'P^' < \T^^t)^^id)^\ >= S+{po,p)p- + S-{po,p)p+ 
-oo 

(2.39) 

S^iP0,P) = i , ^, V .+2 (2-40) 



with ^ ^ 

Tg = T3 COS ^ - n sin — (2.41) 

and T... the time ordering symbol. The state | > in (2.39) denotes the Fermi 
sea of the itinerant holes and the time evolution of V'p is generated by the 
quadratic hole Hamiltonian. 



3. Magnetic Energies with Normal Holes 



In this section we shall calculate the ground state energy and the spin 
wave spectrum at zero temperature using the Hamiltonian (2.11)-(2.17). For 
the sake of simplicity, we assume that only the band et^ is partially filled 
and denote 



with 



^^^+ = ^p = -t'up- Ap- p (3.1) 



= y ^A252 + XStv^cos I + t2^2_ (3 2) 

As we shall see in the next section, the incorporation of possible supercon- 
ductivity orders does not modify the result drastically. We confine ourselves 
within the normal phase of the carriers in this section. 

The basic ingredients of the Feynman diagrams are displayed in Fig. 3 
A dashed line with a plus sign or minus sign represents a magnon prop- 
agator with C = ±1, whose expression is 

i z (3-3) 

for C = ±1. A two point vertex with two incoming or two outgoing magnons 
is associated with ^2iJSv^siv? ^ for C = ±1. A solid line represents a 
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50 100 150 

hole energy (arbitrary unit) 




100 150 
hole energy (arbitrary unit) 



(b) 





(c) 

Figure 2 The density of states of itinerant holes at several canting angles, (a) 6 = 0; 
(b) 9 = 90°; (c) e = 180°. The solid line represents the band with C = 1. and the 
dashed line represents the band with C = — 1. Both kinds of bands degenerate for the pure 
antiferromagnetic order (c) 
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Figure 3 Some elements of the Feynman diagrams of the Hamiltonian (2.11)-(2.17). 

hole propagator and is associated with the expression (2.40). A vertex with 
an incoming and an outgoing holes and an incoming(outgoing) magnon of 

C = 1 is associated with the factor i\\ i^t-{t+) and that with an incom- 



ing(outgoing) magnon of C = — 1 with iAy ^pit_(t+). A vertex with two 
hole lines and two magnon lines of the same C is associated with a factor 
^W^3' ^ Coulomb vertex with two incoming and two outgoing electrons is 

associated with —ijj{D^^^ +D^?^' pi) with q the momentum being transferred 
in the scattering, where the exchange term has been neglected. Other ver- 
tices of (2.11)-(2.17) may also be read off easily but the above ingredients 
are all we need for the following calculations. 

3.1 The ground state energy 

The free energy of the ground state per site in the grand canonical en- 
semble comes from the magnetic terms, (2.12) and the Fermi sea and reads 

/ = 2S^{J' + 2 J cos 0) + ;^ E (3-4) 

p 

Diagrammatically, the Fermi sea term corresponds to ix a closed hole loop. 
The canting angle is determined by the minimization condition 



(I) 



(3.5) 
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and we obtain 

-|-3^E^_''(-.-)- (3.6) 

Diagrammatically, this condition eliminates the sum of tadpole diagrams to 
the order of tree level of the magnon lines and one loop level of hole lines. 
The hole density the number of holes per site, is related to the chemical 
potential through 

"*=(|)e-^E«(-.-)- (3-7) 

P 

The solution of (3.6) and (3.7) determines the canting angle and the chemical 
potential as a function of the hole density. By a Legendre transformation, 
we obtain the canonical ground state energy 

u = f + finh = 2{J' + 2 J cos Q)S^ + JfYl ^i-^p)i^p + ^), (3-8) 

p 



in which the chemical potential and the canting angle are expressed in terms 
of n/j with the aid of (3.6) and (3.7). With a strong double exchange, XS » 
t,t' 

-XS — t'uft — tvftcos — 
2 P P 2 



= — XS — t'u^ — top-cos — — A* (3-9) 



and the sum over p is restricted near the bottom of the band at low hole 
densities n/j << 1, and an explicit expression of follows from (3.6), i.e. 

^°'2 = 4j:s^- ^'-'^^ 

This can be recouncilcd with the formula (1.8) by setting h — riht. Substi- 
tuting this into (3.8) we find 

u = 2{J' - 2J)S^ - {\XS + 2t')nh - (3.11) 

The expression of (3.10) and (3.11) was also obtained by de Gennes [5] with 
the magnetic Hamiltonian (1.3). 
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3.2 The magnon spectrum 



Following the standard strategy, we introduce two self energy functions 
n±(g'o, q) and A-|-(go, q)i where n±(go, Q) represents the sum of all amputated 
one particle irreducible diagrams with an incoming and an outgoing magnons 
of C = ±1 and K±{qQ, q) the sum of all one particle irreducible diagrams with 
two incoming (two outgoing) magnons of C = ±1. Then the full propagator 
with an incoming and an outgoing magnons reads 

-F±{qQ, q)F±{-qo, -q) + A|(go, q) 
and that with two incoming (two outgoing) magnons reads 

V±(*. 5) = -i „, ^ . (3.13) 

-F±{Qo, Q)F±{-qo, -q) + Ai(go, O) 

where 

F±{qo, q) = qo- i^f - U±{qo, q). (3.14) 
The magnon spectrum is determined by the pole condition: 

-F±{qo, q)F±{-qo, -q) + Al{qo, q) = 0. (3.15) 

The derivation of (3.12) and (3.13) follows the same steps as for dressed 
propagators of an interacting boson system with a Bose condensate [9]. 

The self energy functions have been computed to the order of one hole 
loop and the relevant diagrams are exhibited in Fig. 4. 

Note that the Debye screening of the long range Coulomb interaction has 
been included in the C = 1 branch since it plays a crucial rule to stablize 
the canting order. The computation is straightforward and we list in the 
following the main results in the approximation XS » t,t' 

C = 1 Branch[10]: 

U+{qQ,q) = u'^^\qo,q)+U^+\qQ,q), (3.16) 

where n^'*(go,^ denotes (a) + (b), the part without Coulomb vertices, and 
is given by 

A 

2Af 



nf(go,g) = ^E^(-^p)x 



X 



p 
XS 



AS' + go + t'{up - Up-q) +t{vp+ vp^^) cos ® 
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and n^'*(go, Q) comes from the screened Coulomb force, 

ni'H^cg) = (c) = ^c2(go,g)i^(go,g)sin2|, (3.I8) 

where C(^qo,q) corresponds to the couphng between a magnon of C = 1 to 
the density fluctuation and reads 

and 'D{qQ,q) corresponds to the dressed Coulomb propagator 

Diqo,q) = {k)= 3^ (3.20) 
with the density-density correlation function, (T{qo,q), given by 

^(?o, 9) = E — (3.21) 

Likewise, we have 

A+(90, q) = 2JSv^siT? I + Af (50, 9) + A?^ (90, 9), (3.22) 

where ^ 

Af (go, g) = (d) + (e) = sin2 ^ ^ Vf Vf 

V 

X ^ (3.23) 

e^.g— g — On 



and 



Aj^(go, g) = (f) = ^C{qo, q)Diqo, q)C{-qo, -q) sin^ |. (3.24) 
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C = —1 Branch: 



A 



n-(go,g) = ^$]^(-ep) 



XS + qo + t'{u^ - up_^) + t{vp - vp_^) cos § 



2NS 



• 2© V- 
Sin — > 

2 ^ 



-90 + - U^^g) - t{v^_l + V^^^) COS f 



p-2 P+r 

P 2 



go + t'iu^^g - u^_g) - t{v^_^g + v^_q) cos f 



(3.25) 



and 



p+2 p-r 



e 



p+5 P-r 



A_(go, g) = -2JSv^sui' - - 



e 



sm 



2 ^ P+l P-| 

p 



X 



^ P+r 



-go + t'iu^ f — tt^, f ) — t(v^ g- + v^, f ) cos § 
^ p-f p+r p-5 p+r 



go + ^'(Wp.+| - - ^(^^^+9 + I 



(3.26) 



p+5 P-r 

It follows from the equilibrium condition (3.6) for the canting angle that 



F±(0,pz)+A±(0,p£) = 0. 



(3.27) 



The case with p = of this equation is the consequence of the Goldstone 
theorem and (3.27) is found to hold without the large XS approximation. 
Furthermore we find 

F_(0,0) - A_(0,0) = 0, (3,28) 
F+(0, 0) - A+(0, 0) = -1675 sin^ ^ 



sin2 |(< > _ < >2) J] 5(-ep^), 



(3.29) 



where < ... > denotes the average over the Fermi surface 



(3.30) 



A simple expression of the magnon spectrum can be obtained at low carrier 
density Uh « 1 and with a parabolic approximation of the band. The go 
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dependence of n±(go, <?) and A±(go, <?) may be dropped under the condition 
that qo << vpq with vp the Fermi velocity. The second term of (3.29) may 
also be neglected since it is smaller than the first term by a factor of the 

2 

order of nj^ and eq. (3.27) becomes 



^0 



(3.31) 



with qo replaced by u>± we obtain the following long wavelength formula for 
the magnon energy: 



= 4:a± sin ^ ^ / JS 



4:S{t' + t cos 



(3.32) 



for C = +1 and 



{-J'S + ^t')al + JSa'^ sin2 ^ (3.33) 



e 



45 



for C = —1. a± and a of (3.32) and (3.33) are related to the cartesian 
components of q through 



2 



2 / 2 I 2\ 



and 



a 



iql + ql + 2ql). 



(3.34) 



(3.35) 



The expressions (3.32) and (3.33) are different, though not far, from that 
obtain by de Gennes using the magnetic Hamiltonian (1.3). Both branches 

3 

contribute to the specific heat a power law T2 dependence at low T. This is 
different from the observed suppression of the low temperature specific heat 
in the doped samples. It is not likely ,as we shall see in the next section that 
the coupling to a genuine superconducting order will improve the situation. 
An possible explanation is the magnetic anisotropicity, which freezes the spin 
frame to some preferred orientation. The detail of this effect will be discussed 
in the Appendix A. 

The approximate equation (3.31) inspires an effective Hamiltonian for 
low momentum magnon spectrum 
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and ^ 
with 

x± = z.±+n±(0,g-) (3.37) 

and 

y± = A±(0,g). (3.38) 

The magnon of C = — 1 is stable against decay to this order and that of C = 1 
is not. The decay width of the latter is calculated in the Appendix C and is 
found much smaller than the energy (3.32), for large S and small n/j. 

The most remarkable lesson we learned from above Green function for- 
malism is that the role played by the long range Coulomb interaction. If 

the Il^^\qo,q) term of (3.16) were not there, n-|-(q'o,^ would be dominated 

by the second term of 11^ (90) O): which is of the order and is negative. 
This would make (3.29) positive for C = 1 and the corresponding solution for 
magnon energy imaginary. This is associated with the negative coefficient 
of the term of the bulk energy (3.11), which would be unstable against 
clustering if the holes were electrically neutral. It is the long range Coulomb 
interaction that turns (3.29) negative and a real magnon spectrum follows. 
This effect can not be obtained from the magnetic Hamiltonian (1.3). 

4. Monte Carlo Simulation of the effective Magnetic Hamiltonian 

The calculations of the last two sections are restricted to the zero tem- 
perature. For finite temperatures, numerical simulation is the only means of 
exploring the phase transitions. The Hamiltonian (1.1), involving both spin 
and fermion variables is prohibitively complicated to perform a Monte Carlo 
simulation on a sufficiently large lattice. Therefore we switch to the effective 
magnetic Hamiltonian (1.3) with classical spins. The argument which justify 
such an approximation can be found in [6]. 

The partition function corresponding to the Hamiltonian (1.3) is 

Z = J Y[(fSje~^ (4.2) 
j 

where wc have made classical approximation of the spins and the integration 
d?Sj is over the solid angle of Sj with its magnitude fixed. The internal 
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energy per magnetic site is given by 



.2 d 



u^ksT'— In Z (4.3) 
and the specific heat per site is given by 

cv = gf- (4.4) 

The high temperature expansion of the internal energy per site, u, ob- 
tained by expanding the integrand of (5.2) to a power series of and the 
first few terms give rise to 



where 



and 



A = - ^JbS^ + ife^ (4.6) 

A' = J'^S^ - U'b'S^ + lh'\ (4.7) 
5 6 



with G the Catalan's constant given by 

(-1) 



G = y"-^ — ^,. = 0.9160. (4.9) 

n=0 ^ ' 



The corresponding specific heat reads 

-2 ^' + 2^ 16 B 



cv = 



+ 0{{kBT)-^) kB. (4.10) 



.3 (fcBT)2 3 (kBT)^ 

At low temperature, equal partition theorem implies that 

u = UQ + kBT (4.11) 

and the specific heat approaches ks- and therefore can not be trusted there. 

Extensive Monte carlo simulation has been performed on this system. 
There are four parameters in the magnetic Hamiltonian and three dimen- 
sionless ratio can be formed: J' / J < 1, which measure the anisotropicity 
of the antiferromagnetic exchange; b' /h, which measure the anisotropicity of 
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the double exchange, and h/ JS'^^ which measure the strength of the double 
exchange relative to the AF exchange. All these ratios are expected to be 
doping dependent. It follows from the specific heat data for x = that 
J' ~ J, The specific heat data also indicate that significant anisotropicity is 
developed with doping. Therefore, we simulated a number of J' / J ratio for 
b = b' = and only two ratios J' / J = —0.4, 0.5 are considered for 6, b' ^ 0. 
In the latter case, we are more interested in the interplay between the double 
exchange and AF exchange for various b and b' is chosen to equal to b for the 
sake of simplicity. Most of the results except the finite size scaling behavior 
are extracted from a lattice of 2000 sites. The heat bath algorithm was used 
for the simulations without double exchange and Metropolis [11] algorithm 
was used for the simulations with double exchange. 

Fig. 5 shows the phase diagram for J' / J = —0.4,0.5, obtained from the 
MC study as compared with the corresponding mean field result of de Gennes 
[5] . The magnitude of the transition temperature for various b is reduced from 
MF by about a factor 2 but the demarkation point between AF dominance 
and FM dominance remains close to the MF value b/JS'^ = 2.5. The specific 
heat extracted from long MC runs are plotted in Fig. 6 for several b value at 
J' /J = -0.4. 

The two peaks corresponds to the transitions from disordered phase to 
the FM(AF) phase at higher T and the transitions from FM phase to the 
canting phase at lower T. 24000 configurations are included in the average 
for the data points. The estimated error ranges from one percent outside 
the critical region to ten percent inside the the region. As the two transition 
approach each other, its always the peak at higher T that is suppressed. This 
is consistent with the experimental fact that the higher T peak is somehow 
less sharper than the other one [2]. The scaling behavior as the lattice size 
increases from A/" = 250 to A/" = 6750 is displayed in Fig. 7 and the growth 
of these peaks indicate two genuine phase transition as A/" ^ oo. 

The canting order can be viewed as a superposition of FM and AF orders. 
As the temperature is lowered for b/JS"^ < 2.5, the system enter this phase 
from AF order. In another word, a FM order is triggered at the lower 
transition temperature. On the other hand, as the temperature is lowered 
for b/ JS'^ > 2.5, the system enter the canting phase from FM order and a AF 
order is triggered at the lower transition temperature. Since a homogeneous 
magnetic field couple directly to the FM order only, we expect that the peak 
at lower T is sensitive to the external magnetic field for b/JS"^ < 2.5 and 
the peak at higher T is sensitive for b/ JS"^ > 2.5. This is indeed the case for 
the MC results of the specific heat with an external magnetic field shown in 
Fig.8. 
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Figure 5 The Monte Carlo results of the magnetic phase diagrams for (a) ^ = —0.4; 
(b) ^ = 0.5. The solid lines is given by the molecular field approximation [5]. The dia- 
monds denote the onset of an antiferromagnctic order and the crosses denote the onset of a 
ferromagnetic order, x = y = with ks the Boltzmann constant. 

5. Coupling to the Superconducting Order 

Under the assumption that only the band with C — 1 and energy (3.1), 
referred to as the conducting band, is partially filled, the superconducting 

21 




2 



1.8 
1.6 
1.4 
1.2 



A X 



0.8 
0.6 
0.4 
0.2 



0.5 



1.5 



2.5 




(b) 

Figure 6 The Monte Carlo results for the specific heat for several strengths of the double 
exchange and ^ = —0.4; x = and ?/ = c^a in arbitrary unit, (a) = 3 (cross), 



2.75(triangle); (b) 



2(cross), 2.25(triangle). 



pairing, if there is any, has to be of odd parity. The spin orientation of 
the two partners of a pairing with momenta p and -p are parallel, like the 
pairing in the A phase of ii/'e^[12]. But, the orientation itself depends on 
the magnitude of p. The collective mode (Bogoliubov mode) originated from 
the long range order has the C number one and is liable to couple with the 
magnon of the same channel. 
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Figure 7 The scaling behavior of the specific heat for J\f = 250(star), 432(cross), 
2000(square) and 6750(triangle) near the transitions, where ^ = —0.4, = 2.75, x = 
and y = cy in arbitrary unit. 

Introduce the annihilation and creation operators of the itinerant holes, 
atg and with t and s ladling the quantum numbers of the four bands 
(i = ± for C = ±1), the hole-magnon coupling term of the Hamiltoninan can 
be written as 

^3 = -y^ E i^MsC,<^fs^fs + ^MsCH^fs + •••], (5.1) 

where we show explicitly only the terms preserving the C number of the holes. 



p's 



ps 



(5.2) 



with the 2x1 wave function defined in (2.36). Consider the scatering 
process of two holes with opposite momenta 

Hp) + h{-^ ^ h{^) + h{-0) 

with q — j^—p and |g| << \p\. The second order perturbation in i^s produce 
an energy dependent effective potential mediating the two holes, i.e. 



Ap',pAlp, _^[P(go, g) + V(go, g)] 

-A_p. ^^A;. _ .[P(Po, P) + V(Po, P)] 
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Figure 8 The effect of an external magnetic field on the specific heat peaks for ^ = —0.4 
and = 3(a), 2(b). The triangles correspond to zero external field, -f^ = and the 
crosses to =0.5; x = ^fj- and y = cy in arbitrary unit. 



— * 

with Pq = Po + Pq cind P = p + p , which consists of the direct interaction 
(the first term) and the exchange interaction (the second term). For p and 
p' near the Fermi surface and q^O, the direct term dominate and we have 



e 



UeS. = — sin 



(5.4) 
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where u)^ is given by (3.32) and the superscript has been suppressed and we 
have also substitute the approximation 

/2 6 
-tsin-. (5.5) 

at low carrier density. The effective potential (5.4) is rather similar to the 
effective potential mediated by an acoustic phonon, which is retarded and 
attractive, and is felt by holes in all bands. The large core spin S here plays 
a similar role as the large ion mass for the electron (hole)-phonon coupling 
and we expect that the Migdal theorem [13] holds for sufficiently large S. 
Furthermore such an attraction would not be there for just ferromagnetic 
order alone, = 0, a point stressed also in [14]. There is a simple physical 
reason: Like the electron-phonon coupling, the effective attraction is achieved 
through local polarization of the ions by charge carriers, the attraction here 
is caused by local distortion of the ion spins. The counter reaction of this 
distortion would flip the carriers spin, which would cost huge amount of 
energy ~ A^" with only ferromagnetic order, since the spin orientation in 
the conduction band point to a unique direction. With canting order, the 
spin orientation in the conduction band varies with the Bloch momentum 
and the local ion spins can be distorted at much less cost of energy. In 
addition to (5.4), there may be other attractions induced by large hole- 
phonon coupling caused by Jahn-Taylor effect. On the other hand, we are still 
far from inferring that there are sufficient attractions to beat the Coulomb 
repulsion. Furthermore the exchange term in (5.4) has to be considered for 
single specie pairing. 

Leaving aside the pairing mechanism, we may still study the coupling of 
the magnon with the superconducting order by by adding a pairing term to 
the Haniltonian (2.11). The total Hamiltonian reads now, 

U = H + Y.^\^'i^ (5-6) 

where ttp*, at stands for the annihilation and creation operators of the itiner- 
ant holes, and /Jg', 
Cooper pairs, i.e.. 



ant holes, and (3fj, pt stands for composite boson operators which represent 



f^? =77^T! 9p,qap+f-p+l (5-7) 
p 

with Yl'p extends over half of the BriUouin zone only. The pairing wave 
function, g^^^^, is odd under the space inversion, i.e. 

9p,q = -9~p,q (5.8) 
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on account of the anticommutation relation among aps [15]. The supercon- 
ductivity is implemented through the condensation of the pairing operator 
of zero total momentum, i.e. 

< S\P^=o\S >= y/NB, (5.9) 

where B is the long range order parameter and is chosen to be real and 
positive for the sake of simplicity. At the presence of (5.9), the part of the 
total Hamiltonian (5.6), which is quadratic in a and a) reads 



Hsc = ^ ^pO^jf^p + XI ^^i^p^-P + V-p^^) (^-10) 



p p 

with Cp* given by (3.1) and 

Hsc can be diagonalized by the Bogoliubov transformation 

aff = aff cos Off — e^'^pa^ ^sin 9ft 
p p p —p p 

= e~*'^PQ;p'sin 9^+ q;1-.cos 9^, (5.12) 



where 



cos2^p^=g 

ISA 

sin2%=i^ (5.13) 
op 



with £p= \/^p+ and we have 



Hsc = lYi^p-^p) + Yl ^P^I^P- (^-14) 



2 

p 



The free energy density of the ground state, (3.4), is replaced by 

1 1 
Jf2 



f = 2S\J' + 2JCOS + - ^p)- (5-15) 



The canting angle is determined now by 
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and the expression of the carrier density (3.7) becomes 



In addition to the strong double exchange and low carrier density, we assume 
that the gap energy comparable or smaller than the Fermi energy. It follows 
then that the expression (3.10) for © remains an approximate solution of 
(5.16) and (5.17). The reason is that the canting is induced by the delocal- 
ization of the charge carriers, which lower the energy of the kinetic energy by 
an amount of the order hopping amplitude, while the approximation amounts 
to assume that other energy scales (Fermi energy and gap energy) are much 
smaller than that. 

The superconducting order affect the magnon spectrum in two ways. The 
first of them amounts to dress the hole lines in Fig. 4 with the long range 
order. On writing the C = 1 branch of the hole propagator, S'^(po,p) of 
(2.40) as 

sHpo,p) = + — %iT7r)' (^-is) 

V nn — »n — + 2 A^?/ 



where Zi are 2x2 matrices 



Only the first term is subject to the modification. The factor 



Po - 

of (5.18) ought to be replaced by 



= / dte'P°^ < \Tap{t)a^Oy\ > (5.20) 

J —oo 



r dte'Po* < S\Tajj{t)ap{0)^\S >= (5.21) 

J-oo Po ~ 

where \S > denotes the ground state with the long range order, (5.9), and the 
time evolution of ap{t) is governed by the Hamiltonian (5.14). Anomalous 
propagators 

/ dte'P°^ < S\Tap{t)a_p{0)\S >= ^ ^ (5.22) 

J-oo Po — ^ig 
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and 

dte'P°* < S\Ta_^t)^a_^0)^\S >= ^ (5.23) 

-oo Po — 

should also be included. The hole propagator of the C = —1 channel 
S~{po,p), however, remains unchanged. With these modifications, the same 
magnon spectrum (3.32) and (3.33) emerge as the leading order of the ap- 
proximation stated following (5.17). 

The second type of influence of the superconducting order is the mixing 
of the magnon of C = 1 and the collective (Bogoliubov mode). Instead of an 
extensive diagrammatic analysis, we take a short cut. We first calculate the 
amplitude of the following two processes: 

I: \i >= Cl\S >, 1/ >= pI\S >; 

II: \i >= Cipl^lS >, 1/ >= \S > 

where, \i > and |/ > denote the initial and final states. Using (4.1), 
(4.7), we find that 

« . lim < m,]i >= ^ =in I E (5.24) 

P P 

From now on, we regard /Jg* and /3t as elementary boson operators and 
introduce the following effective Hamiltonian for the mixing 

^eff. = E ^P<i<P + E' VpiCpC-p + Clp4) 
P P 

+ E + E' vpiPpP-p + 
p p 

E (f^P^P + f^P^-P + ) ' (5-25) 
p 

where and here are the same as xt and y^ in (3.38) and (3.39) and the 
structure bosonic part Hamiltonian can be obtained simply from the boson- 
fermion model of superconductivity [16]. The low momentum limit of the 
coefficients x^, y^^, and 7]^ are all positive and satisfy 

xp-yp=sp, (5.26) 

Xp + yp = c^, (5.27) 
^p-Vp = s'p (5.28) 
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and 



_ ^6 



(5.29) 



with Sp* and s'p two positive quadratic form of the momentum p. The posi- 
tivity of Tj^ comes from the repulsive nature of the boson-boson interaction 
and the choice of the condensate phase made in (5.9). The eqs. (5.26) and 
(5.28) are the consequences of the Goldstone theorem and eq. (5.29) foUows 
from the Coulomb interaction with ojq the plasma frequency. The set of 
Heisenberg equations of motion read 



Kp = xpCp + ypC-p + i^{Pp + P-^), 

-iClf = x^lp + VpCp + «(/3p + Plf), 

0p = ^pPp + VpPlj^ + « (Cp + Clp-) , 



CpPlp + VpPp + <Cp + Cl 



-i/3l 

With the ansatz that cc e 
obtain the secular equation for the energy Q, 



-int 



C-f! oc e*"^ (5p oc e-*^« and (3. 



X 



p 
Vp 

K 
K 



Vp 

K 
K 



K 



Vp 



K 
K 

Vp 

^p + ^ 



(5.30) 

(5.31) 
(5.32) 

(5.33) 



OC e*"% we 



(5.34) 



and its solutions: 



y} + fp 



v} 



The stability requires that > which is equivalent to 

{xp + yp){^p + Vp) > «^ 



(5.35) 



(5.36) 



and is amply guaranteed by the relations (5.27)-(5.29) for p 0. It follows 
from (5.26)-(5.29) and (5.35) that the low-lying spectrum is given by 



c'gp + Oip^), 



(5.37) 



where the O(p^) represents the correction from the superconducting order 
and is therefore rather small. 
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6. Concluding Remarks 



We conclude this paper by a few comments related to the real system, 
Sr2Y Rui-.xCuxOq and Ba2Y Rui-xCuxOq. 

The Monte Carlo simulation and spin wave fitting of the specific data 
[2] for X = indicates that J — 1.28 ~ 2.06meV. The closeness of the 
locations of the two specific heat suggests that we are in the neighborhood of 
the merging point of the three phases (FM(AF), canting and disorder) and 
implies -j^ = 2.5, which together with (1.8) gives rise tot = 36.0 ~ 57.9meV. 
The large double exchange approximation is therefore quite good for the 
real system. Assuming only one component of the d-orbit of Ru^'^ gets 
delocalized, there would be two states per site. Making an equal partition 
of these delocalized states into four bands, there would be 0.5 states per site 
and per band. With only one band filled partially, the filling fraction of the 
band would be 4x. There the low density approximation is not a good one 
for the observed samples (with x=0.07, 0.10). With all three d components 
delocalized and degenerate, the filling fraction can be lowered to 4a;/3 and 
the accuracy of the low density approximation can be improved. 

So far we have concentrated our attention to the magnetic properties of 
the double exchange system. Now we come to the superconductivity of the 
system. Since it has not been rule out that observed superconductivity is 
caused by the cuprate impurities, we shall list in the following the pros and 
the cons to a genuine superconductivity of a pure system. 

1) Imagine the system Sr2Y Rui-xCuxOq consists of several domains of 
pure phases of the double perovskite Sr2Y Rui-xnCux^ and cuprate Y Sr2Cus-yf,Ruyf,OjSm-^ 
Suppose that the size of each domain is sufficient large that the interface ef- 
fect on the stoichiometry may be neglected, we expect that 




(6.1) 



n 



k,m 




(6.2) 



n 



k,m 




(6.3) 



n 



k,m 



and so on, where is the total number of Y ion, Nn is the total number 
of Sr2YRui-xnCux„0e molecules and N^.m that of Y Sr2Cu'i-yj^Ruyfij_5,^ 
molecules. The eqs (6.1)-(6.3) is the consequence of the balance the numbers 
of F, Sr, Ru and Cu ions. Combining (6.2) and (6.3) we find N — + 
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•^^km^k,m- It follows from (6.1) then that Nj^^^n = 0) which rules out the 
bulk impurity phase of cuprates. 

2) The coherent length of a cuprate superconductor is comparable to 
the lattice spacing. A grain with a small number of unit cells of cuprate is 
already able to sustain a long range order. It may happen that a large number 
of such grains of cuprates get embedded inside the bulk double perovskite 
sample. On the other hand, in order to support the bulk resistivity transition, 
Josephson tunneling between different grains has to be significant. This 
process, however is made rather difficult by a canting magnetic order since 
the splitting of the spin degeneracy in the bulk material reduces drastically 
the tunneling probability of a singlet Cooper pair. 

3) If the observed superconductivity is a genuine one coexisting with the 
canting magnetic order, the pairing is between the itinerant holes with oppo- 
site momenta and is of odd parity. The doped Cw-ions plays the role of impu- 
rities, in addition to provide itinerant holes. (In this aspect, Sr2YRui-xCuxOQ^ 
differs from the cuprates. The doping of the latter does not alter the Cu02 
planes.) These impurities offset the degeneracy between opposite momen- 
tum and thereby suppress the superconductivity. This is similar to the effect 
of magnetic impurities on the singlet pairing. If this factor dominates, the 
transition temperature would be lowered with doping, opposite to the obser- 
vation. 

Odd parity pairing would be a distinct feature of the superconducting 
order with canting magnetic order and the parity symmetry is spontaneously 
broken, i.e. 

V\S>^\S>, (6.4) 
where V stands for the ordinary inversion operator defined by 





(6.5) 




(6.6) 




(6.7) 



and = 1. Unfortunately, such a parity broken does not produce any 
observational consequence for a homogeneous super phase. A new operator 
which combines the inversion and a charge rotation through |, 

V' = e*5^P (6.6) 

with Q the charge operator, leave the condensate invariant. Acting on the 
observables bilinear in hole operators, momentum, electric current, etc., V' 
is completely equivalent to V since these operators carry no charges. The 
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square of P', however, assumes different eigenvalues for charge even and 
charge odd states, V'^ = {—)^- But the states with Q =even never coher 
with the states with Q — odd, even with a condensate of Cooper pairs of 
charge 2e. 

With both an odd parity superconductor and an even parity supercon- 
ductor coupled together, none of V or V' will leave the ground state invariant 
and observational consequences of the odd parity pairing can be extract. 

Furthermore, the pairing state \S > breaks also the time reversal invari- 
ance and this will lead to additional observational effects which may serve 
as an evidence of a genuine superconductivity. 

Another possible way to distinguish a genuine superconductivity from 
that caused by cuprates is to examine the anisotropicity of the observed 
superconductivity. X-ray data suggested only moderate anisotropicity of 
the crystal structure in contrast to the strong anisotropicity of the cuprates 
because of the quasi-two- dimensional crystal structure. 

In addition to the magnetic properties presented in this paper, we have 
also calculated the usual optical response functions on account of the exotic 
band structure. The motivation came from the A term of the Hamiltonian 
(1.1), which is, with a ferromagnetic order canting order, equivalent to an 
ultra-strong magnetic field acting on the spin degrees of freedom. The typical 
Farady angle associated to the optical activity because of T-violation of the 
magnetic ordering turns out to be about 10~^ radians per wavelength in the 
visible region and is characterized by a different wavelength dependence from 
that of the ordinary magneto-optical effect. The details of this work will be 
reported in a forthcoming publication [17]. 
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Within the orthorhombic symmetry, we assume that z-axis is the easy 
axis of the ferromagnetic order and x-axis is the easy axis of the antifer- 
romagnetic order. Wc add the foUowing anisotropy terms to the effective 
magnon Hamiltonian (3.36) 

-f^mag. = -f^mag. + ^ ^ [^lb+ ~ {^ab+ ' ^)^] + ^ ^ [^b- ~ 0ab- ' ^)^]) 

<ab> <ab> 

(A.l) 

where 

Sab± = l{Sa±Sb) (A2) 

and the summation extends to all inter-layer bonds. The new terms leaves 
C-parity intact. Applying the large spin expansions (2.9) and (2.10) with 
C = z and ^ = X, and making the Fourier transformations (B.l) and (B.2) 
below, we find that 

+ E ^9^'4^? + E' y^(^<i^-? + ^-f4) ' ^^-3) 

9 



where 



X- 



and 



yf = V^- T ^g) sin^ | - ± v^) + ^w;_(4 ^ ^;g-.) sin^ | (^.5) 

with and given in (3.37) and (3.38). Goldstone theorem ceases to hold 
and energy gaps begin to open between the ground state and the excitations. 
The energy gap for the magnon of C = 1 reads 

(5+ = ^^2JSw- sin ^ {A.Q) 
and that for the magnon of C = —1 reads 

5- = "^s^w^ {wj^ cos^ ^ + sin^ . (^.7) 

With a fixed canting angle 0, the spin frame can be viewed as the body frame 
of a rigid rotator and its orientation can be specified by three variables, say 
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Euler angles. On the other hand, for smaU vibrations about the laboratory 
frame (i.e. lattice frame here), the dependence on the precession angle is of 
higher order in magnitude. This is why only two parameters w± are included 
in (A.l). 



Appendix B 



1. Large spin expansion 

In this part of the appendix, we shall display the steps which lead from 
(1.1) to the expressions (2.11)-(2.17). Making Fourier expansion with respect 
to the lattice Bloch momentum p, 



P 



(a) ip-Ra 

P 



P 



{B.2) 
(S.3) 
{BA) 



with a labeling the sites on A and b labeling the sites on B. Substituting 
(2.9)-(2.10) and (B.1)-(B.4) into the Hamiltonian (1.1), we obtain 



where 



H = Hq + Hi + H2 + Hs + H4 + Ucoul.i 

Ho = 2A/'S'2(J' + 2Jcose), 
Hi = AJS^VH{cq + 4 - do - 4) sin © 
H2 = 2J'Sj2i-'^ + Up)i4^P + dyip) 



{B.5) 

{B.6) 
{B.7) 



+2JS 



e 



- 4 ^(ctc^ + did^) cos B + ^ v^{cidp + dlc^) cos^ — 
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+2 J \ sin Q'S^ Vff\ — Cffd^ ^d^ q — ctd^ 



P+i 



7^4 



2j; 



' P+I P-l 1 p-i P+IJ 
p'-f P+I f+i P-l f-l P+I P'+I P-|/ 



(S.9) 



p,p',g 



P,P ,9 



p+§ P-J 



A 



^-^^ \ F-| P +5 P+I P 5 P+^ P 5/ 



p-i f-i P+i p+r 



(5.10) 



2 

jVejT^^L « p+l p'-l P +5 P~2 



P,P >9 



In the above expansion, we have combined hole operators ^^^^ and ^^^\ each 
of two components, into a 4 x 1 column matrix. 



*p 



(a) 
P 
(6) 



(S.12) 



Accordingly, r's and p's stand for the 4x4 Dirac matrices (2.18)-(2.19) with 
T = Ti^ + T2'ri + TsC and p± = The 4x4 matrix E^^ is given by 

= -t'ujj - tv^pi + ^XS{p+Ca ■ r + P-Cb ■ r) {B.13) 
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and the unit vectors e^^.^^ is defined as 



1 

and 



tip- and vp in the above equations are given by (2.25) and (2.26) 

In addition to invariance under the tetragonal crystal rotation, inversion, 
the following transformation C, 

CcpC"^ = -dp, {B.IQ) 



Cdf,-^ = -cp (B.17) 

and 

C^pC-^ = C^p (S.18) 

with C — —ipie~^'^^ = — Pits also leaves the Hamiltonian (B.5)-(B.ll) 
unchanged. This transformation amounts to interchange the sublattices A 
and B while making a tt rotation of the ion spins about 2;-axis. The time 
reversal invariance is broken by the magnetic ordering. This symmetry can 
be diagonalized via the following unitary transformation: 



<'p = ^i^P + Cp), {B.19) 



where 



^p=Vijp, (5.21) 

F = (p_ + p+Ts) (pi + p3) ( COS ® - iT2 sin ^) . {B.22) 

In terms these new operators, the transformation laws (2.29)-(2.31) follow 
and a simpler form of the Hamiltonian, (2.11)-(2.17), emerges. 

B.2, Coulomb vertex 

In what follows, we shall calculate the Coulomb vertex. Labeling the 
sites on the sublattice a and those on the sublattice by b, we have 

r 

Ucoul. = ^ \U (^ : QaQa '■ + ^ ' QbQb ■) 
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: QhQi/ : 



where 



2E 



a6 



: QaQb ■ 
Ra — Rb\ 



(5.23) 



(S.24) 



with j = a,b and :(...)• denotes the normal ordering. Making the Fourier 
transformation in (B.3)-(B.4), we obtain 



Ucoul. 
where 

and 



2Afe 



^ ^ P+§ P-i ^ 



(S.26) 



(S.27) 



4-= E*^-^?^3*--«~ = E^^-^?Pi^zr-« 

with and V'p' related by the transformation (B.21). and S'^ denote two 
Madelung type sums, 



1 Q-iq-Rn 



n/0 



and 



with 



and 



1 x-^ e~*^'-^" 



3 47r 



I 



Rn = ^[ni{x + y)+ n2{-x + y)] + n^lz 



(5.28) 



(S.29) 



(S.30) 



[ni + ^) (:r + y) + (na + ^) {-x + y)] + (rig + ]^lz (S.31). 



It follows that gt is C even and is C odd. The whole Coulomb energy is 
therefore invariant under C, but both branches of the magnons couple to the 
Coulomb interaction. and D^2'^ in the text is given by 

Df = S^+S'. (S.32) 
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and 



(S.33) 



The calculations of D^S^ and D^S^' at low q are similar and we illustrate 

the method by D^^^ . We first exponentiate the denominator of the summand 
of (B.28). 

1 _ V2 

Rn I 

and interchange the order of the integration and the summation, i.e. 



Jo 



where 



U + 



V2 
47tI 



Ct2 



as 



2it 



d2,{z\T) = I + e*'''^^ cos 2nz 
is one of Jacobian theta function [18]. Using the formula 



(S.34) 

(S.35) 
(S.36) 

(5.37) 



we have 



3 ^"±+"3 , / ai,t 

+t-2e STF^-i^g 



-i 



a2 1 1 
2t 't 



V2 



-i 



4t^2t 



{B.38) 



Only the "1" term of each -ds factor of (B.38) is sensitive to the limit ^ — > 
upon integration. By isolating out this term and setting ^ = for the rest, 
we arrive at 



U + 



where 
C-- 



-f 



t~ 



4 



AttI 

Likewise we obtain that 



l(2a\ + al) 



+ C, 
21 



0.27214 



l{2a\ + al) 



0.092711 



(S.39) 



(S.40) 



(S.41) 



Substituting (B.39) and (B.41) into (B.32) and (B.33), we end up with (2.27) 
and (2.28). Only the magnon of C = 1 couple to the long range tail of the 
Coulomb force. 
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Appendix C 



In this appendix, we calculate the decay rate of the magnon with C = 
1. Neglecting possible couplings to the superconducting order, the effective 
magnon Hamiltonian for C = 1 corresponds to the first two terms of () with 



and 



a;p^=iy+ + n+(0,p) 



= A+(0,p). 



(C.l) 



(C.2) 



A mass-shell magnon operator is given by the following Bogoliubov transfor- 
mation: 

7p = Cp cosh Xp + Cl^sinh Xp, 



' — r 



Cp sinh Xp + Cl^cosh Xp 



(C.3) 



where 



coshxp 



X 



sinhx^ 



+ 



with ujp = y/x?,— whose low momentum limit is given cxplicity in (3.32). 

The amplitude for a magnon with momentum q to decay into a pair of 
an electron and a hole with momenta p ± f near the Fermi surface, 
corresponds to the superposition of the two sets of diagrams in Fig. 9, i. e. 

-^M = ^^li ((Fig.9a) cosh Xq - (Fig.9b) sinh X?) (C.4) 

With the aid of the Feynman rules listed in the section 3, we find that 



a- cosh Xq - -7+ sinh Xq 



+2- 



(0) 



^ sin ^[C(a;,g) cosh Xg- C{-lo, -g)sinhxg] 



p-i+ 



(C.5) 

where the functions D^^\ a{uj,q) and C{uj,q) are defined in section 3 with 

LJ = Ujt. 



39 



Following the same reasoning employed in section 3, we approximate 
C{u>, q) and (t{u>, q) by C(0, q) and a{0, q) and obtain, after some algebra 



2AfS 



-\-{vp- < V >)2(cosliXg - sinhxg-)^ sin^ ^ (C.6) 

for small q, where < .. > denotes the same Fermi surface average defined in 
(3.30). The decay width is given by 



With parabolic approximation of and and low carrier density, n/j << 1, 

2 

is of the order of so is the Fermi energy measured from the bottom of 

the conducting band, ep = + 2t' + 4tcos ^. The order of the first term of 

2 

(C.6) is nj^, the order of the second term is n/j and that of the third term is 

4 

nj^. Keeping only the first term of and completing the phase space 

integral for small q, we obtain 

2 J Stt''^e pa sm^ (5 cos^ (5 6.06 , 

Vff = — , tan — sm — , (6.8) 



{f + 1 cos f )2 (t' sin^ /5 + t cos I ) i y^tcos 



2 



where we have parametrized the magnon energy ujt by ojt = ca sin /? with 

a = ■sj af + a2 + a"^ and sin/5 = a±/a. The contribution from the second 

term in the bracket of (C.6) has been neglected in (C.8) since it is of higher 
power of the hole density, n/j, when n/j << 1. As far as the order of mag- 
nitude concerned, we have c ~ JS, t' ~ t, and JS^ ~ tn. It follows then 
that 

J>in|«l (a9) 
for large spin and small carrier density. 
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(b) 

Figure 9 The Feynman diagrams for the decay of a magnon with C = 1. The open circle 
with four legs stands for the dressed Coulomb vertex as in Fig. 4. 

REFERENCES 



1. M. K. Wu, et. al., Z. Phys., 102, 37 (1997). 

2. M. K. Wu, et. al., submitted to Phys. Rev. Lett.. 

3. C. Zencr, Phys. Rev., 82, 403 (1951). 

4. P. W. Anderson and H. Hasegawa, Phys. Rev., 100, 675 (1955). 

5. P. G. de Gennes, Phys. Rev. 118, 141 (1960). 

6. A. J. Millis, P. B. Littlewood and B. I. Shraiman, Phys. Rev. Lett., 74, 5144 
(1995). 

7. P. D. Battle and W. L J. Macklin, J. Solid State Chem., 52, 138 (1984). 

8. T. Holstcin and H. Primakoff, Phys. Rev. 58, 1098 (1940). 

9. N. M. Hugenholtz and D. Pines, Phys. Rev. 116, 489 (1959). 

10. M. Creutz, 'Quarks, Gluons and Lattices', Cambridge University Press, 1983. 

11. To avoid cumbersome expressions, we have suppressed the infinitesimal 
imaginary part of the denominators of 11^^ {qo,Q), {qo,Q), C{qo,q) and 
a{qo,q) and supplemented with an equivalent rule that the imaginary parts of 
these functions for real go has to be negative, as is evident from Dyson- Wick 
perturbation theory. 

12. A. J. Leggett, Rev. Mod. Phys., 43, 331 (1975). 

13. See for example, A. L. Fetter and J. D. Walecka, 'Quantum Theory of Many 
Particle System', Chapter 12, (McGrawHill, New York, 1971). 



41 



14. W. E. Pickett, Phys. Rev. Lett., 77, 3185 (1996). 

15. Rigorously speaking, pairing term in the Hamiltonian (5.6) violates the 
tetragonal symmetry of the original Hamiltonian (1.1). On the other hand, 
the Goldstone theorem can be verified to the leading order of the low carrier 
density. A rigorous treatment of the pairing effect should respect the original 
symmetry and therefore should be extended to all bands. 

16. R. Friedberg, T. D. Lee and H. C. Ren, Phys. Rev. B42, 4122 (1990). 

17. H. C. Ren and M. K. Wu, in preparation. 

18. E. T. Whittaker and G. N. Watson, 'A Course of Modern Analysis', 4th ed.. 
Chapter XXI, Cambridge University Press, 1927. 



42 



